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Abstract: When ultrafast noncritical cascaded second-harmonic genera- 
tion of energetic femtosecond pulses occur in a bulk lithium niobate crystal 
optical Cherenkov waves are formed in the near- to mid-IR. Numerical 
simulations show that the few-cycle solitons radiate Cherenkov (dispersive) 
waves in the A = 2.2 — 4.5 /xm range when pumping at Ai = 1.2—1.8 /im. 
The exact phase-matching point depends on the soliton wavelength, and 
we show that a simple longpass filter can separate the Cherenkov waves 
from the solitons. The Cherenkov waves are born few-cycle with an 
excellent Gaussian pulse shape, and the conversion efficiency is up to 25%. 
Thus, optical Cherenkov waves formed with cascaded nonlinearities could 
become an efficient source of energetic near- to mid-IR few-cycle pulses. 
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1. Introduction 



Nonlinear optics dawned when second-harmonic generation (SHG) was demonstrated 50 years 
ago H]. Soon after, it was predicted Q and experimentally demonstrated IS that cascaded SHG 



could lead to a Kerr-like nonlinear action. The cascading occurs when the harmonic generation 
over a propagation length L is strongly phase-mismatched |AA:L| ^ 1: only a fraction of the 
fundamental wave (FW) is up-converted to the second harmonic (SH) after a coherence length 
;r/|M|, and after another coherence length it is back-converted to the FW. Due to the different 
FW and SH phase velocities when Ak^O, the back-converted FW is phase shifted relative to the 
unconverted FW. On continued propagation the FW effectively experiences this cascade of up- 
and down-conversions as a Kerr-like nonlinear refractive index change ziw = "casc^ proportional 
to the intensity /. Since nj^^^ °^ —d^ff/Ak the phase-mismatch Ak = k2 — 2k\ controls both the 
magnitude and sign of the cascaded nonUnearity This has been used to study, e.g., spatial, 
temporal and spatiotemporal solitons |l4l|5]|6l, high-energy pulse compression fTl [8l l9l [TOl . 
supercontinuum generation ITTl . and all-optical signal processing ||T2l . 

Recently we investigated high-energy pulse compression through few-cycle solitons gener- 
ated in cascaded SHG, and showed that the soliton can couple to dispersive waves ifTSl [T4l . 
Dispersive waves were first predicted as a consequence of perturbing a stable temporal soliton 
in the nonlinear Schrodinger equation (NLSE) fTSll . The first experimental observations came 
shortly after in a mode-locked dye laser lfT6l and in single-mode fibers ifTTlfTSl : both systems 
are well-described by the NLSE. The dispersive wave is linear in nature and Ref. |191 showed 
that its spectral location is a result of a phase-matching condition to the soliton, and that its 
strength relates to the overlap to the soliton spectrum. There it was also pointed out that the 
dispersive wave shed by the soliton is reminiscent of optical Cherenkov radiation, and it has 
recently found applications in supercontinuum generation 1201 . for ultra-short pulse synthesis 
II2TI . broadband frequency combs [22] and UV femtosecond pulse generation Il23ll24ll . 

Cascaded SHG can induce both self- and cross-phase modulation cubic nonlinear terms 1251 
|26l |27l . and in the limit of a weak SH the FW is well-described by an NLSE, just like the 
cases above. In Ref. lT4l we used this property to show that in self-defocusing cascaded SHG 
the Cherenkov wave is emitted at longer wavelengths than the soliton. Self-defocusing solitons 
namely form in the normal group-velocity dispersion (GVD) regime HI, and the Cherenkov 
wave is emitted in the non-solitonic anomalous GVD regime at longer wavelengths. 

Important C-H, N-H and O-H stretching modes reside around A = 3 /im, and probing or 
controlling these require access to energetic few-cycle mid-IR pulses. Unfortunately, the beta 
barium borate crystal studied in Ref. lT4l has an absorption edge around A = 3 /xm. Recently we 
showed experimentally that noncritical cascaded SHG in lithium niobate (LN) can support near- 
IR few-cycle solitons ITOl . and the predicted phase-matching to mid-IR Cherenkov waves lfT4l 
was confirmed numerically. Here we give a detailed numerical investigation of these mid-IR 
Cherenkov waves: they lie in the A = 2.2 — 4.5 /im range (the exact resonance point depends on 
the soliton wavelength and strength), they can be long-pass filtered, have excellent pulse quality, 
and are born with almost transform-limited few-cycle duration. The conversion efficiency is 1- 
25%, depending on how far the Cherenkov wavelength is from the soliton wavelength, implying 
that this scheme could be used as an efficient source of energetic few-cycle mid-IR pulses. 

2. Theory 

In Ref. lfT4l we used that cascaded SHG with an effective self-defocusing nonlinearity under 
suitable conditions can be reduced to a normalized NLSE for the FW field Ui 
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c c 



The effective soliton order A^gff 1281 appears due to the following normalization choices 
^ = z/Ld^i, t' = r/Fin, where is the input pulse duration, Ld,i = ^in/l^f is the dispersion 
length, and is normalized to the input intensity /in. The normalized dispersion operator 



is b\ = l^2=2i"'Sr7^, where 5^' = Lo.ikf {T^m\)-\ and kf = d'"ki{co)/d(o"'\co=co,. 
Interaction in a bulk medium implies ki{co) = ni{co)co/c, where ni{(o) is the FW linear re- 
fractive index. The cascaded nonlinear response is n^.^^^ = —2(Oidlff/[c^eon\n2^k] 1141 . and 
"iCerrei electronic Kerr nonlinear refractive index. Besides being in the cascading limit 

(|AA:L| » 1) with a self-defocusing cascaded nonlinearity (Ak > 0), the conditions for using the 
above model to describe the cascaded interaction are |fT4ll (a) a broad-band cascaded nonlinear 
response, which requires Ak > ^fj2/2^2 ^^^^ so-called stationary regime) 1291 . (b) a modest 
effective soliton order, which serves to reduce self-steepening effects (either direct or through 
cascading |l3Ql), (c) a weak SH, in order to minimize cascaded IZSl l and Kerr [|28| cross-phase 
modulation terms, and (d) a negligible Raman response (for the LN crystal we consider in this 
paper we calculated the characteristic Raman time Tr = 0.4 fs, which justifies this assumption). 

We seek an overall effective self-defocusing nonlinearity n'^g < so solitons can be ex- 
cited in the normal GVD regime. In the simplest case a soliton will shed Cherenkov radiation 
according to the phase-matching condition ^dw(ft'dw) ~ ^soi(t'^dw) ~ 0. In cascaded SHG the 
FW forms a soliton at some frequency Osoi- The dispersion relation reflects its nondispersive 

nature: ^soi(«) = ^i(«soi) + (« - COsoi)/vg^soi + ^soi, where v^.soi = l/^l''(«soi) is the soUton 
group velocity and q^oi is the soUton wave number. The Cherenkov wave dispersion is simply 
determined by the FW wavenumber kdw{(o) ~ {(o)- This gives the phase matching condition 

II20I ki (cOdw) — ki (©sol) — (cOdw — 0}so\)k[^\(Oso\) — qso] = 0. The soliton wavenumber can be 
estimated from Eq. ([TJ as q^oi = «iff^soifi'soi/2c fl^. Predicting accurately the soliton intensity 
/sol is difficult, but for low soliton orders its contribution is minimal (see Ref. JM)). Therefore 
the phase-matching condition can be approximated by ifTSl 

^l(«dw) -fcl(«sol) - («dw - COso\)k[^\cOso\) = (2) 

In Fig. [Ha) the phase-matching condition Eq. ^ is shown for 5% MgO doped LN (MgO:LN) 
X-cut (9 = n/2) for noncritical (type 0, ee — > e) interaction. Pump wavelengths of Ai = 
1 .2 — 1 .8 /im should generate Cherenkov waves in the desired region around A = 3 /im. A pre- 
requisite is exciting a soliton broadband enough to spectrally overlap the Cherenkov resonant 
wavelength, and in Ref. (TOl we demonstrated this experimentally: a few-cycle self-defocusing 
soliton was formed in LN through noncritical cascaded SHG pumped with 50 fs pulses at 
Ai = L3 /im. The key factor is a large and ultrafast self-defocusing cascaded nonlinearity and 
we showed that it remains so at other pump wavelengths as well. Consider Ai = 1 .65 jUm (the 
case shown below): despite a huge phase mismatch, Afc = A:2 — 2ki = 283 mm^' corresponding 
to 11 /xm coherence length, the strong effective nonlinearity t/gff — 20 pm/V ||3T1 is able to 
create a large negative cascaded nonlinearity n'^i^Q — —40 x 10^^" m^/W without using quasi- 
phase matching. A large group-velocity mismatch (GVM) c/12 = —250 fs / mm at Ai = 1 .65 jUm 
gives a just 120 /im walk-off length for a 50 fs pulse, which should render femtosecond cas- 
cading inefficient, but the large phase mismatch ensures an ultrafast cascaded nonlinearity: it is 
of nonresonant nature (practically instantaneous with a sub-fs response time) as the interaction 
occurs in the so-called stationary regime ||29l Ak > Ak^r = c/j2/2^2^' — 90 mm^'. 

3. Numerical simulations 

The plane-wave numerical simulations below are based on the slowly evolving wave ap- 
proximation II30I l28l . Since there are no direct measurements in the hterature of the elec- 
tronic and Raman Kerr nonlinear strengths of LN (see discussion in Refs. 1321 [TOl ). we use 
"Kerr = 30 X lO^^O ^2/-^^ ^t A 1.30 ^m, and a Raman fraction /r = 52% HU). Thus, 
"kerr el = ~/s)''kerr = ^"^'"^ ^ 10^^" m^/W, making the effective nonlinearity self-defocusing 
for Ai > 0.8 /im ifTol and solitons can form in the normal GVD regime (A < Azd = 1 -92 /im). 
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Fig. 1. (a) Cherenkov phase-matching condition vs. soliton wavelength, Eq. for cas- 
caded noncritical SHG in MgO:LN (using the Sellmeier equation at 24.5°C from Ref. 1331 ). 
Data from Ref. 1141 . (b-d) Simulation of noncritical cascaded SHG in a 20 mm MgO:LN 
crystal using Ai = 1.65 /im, 50 fs FWHM, A'eff = 3.0, /i,, = 77 GW/cm^. (b) Intensity and 
pulse duration after longpass filtering the FW (A > 2.5 jlm transmitted), (c) The filtered 
Cherenkov wave at various propagation stages (the conversion efficiency is indicated above 
each stage), (d) shows the FW time trace and spectrum. In the spectrum 'A' and 'N' mark 
anomalous and normal dispersion regimes, and the color scale is saturated for clarity, (e) 
XFROG spectrogram (linear scale) at z = 20 mm calculated with a 10 fs gating pulse. 



Figure [U (b-d) shows a numerical simulation using Ai = 1.65 jUm 50 fs FWHM Gaussian 
input pulses with A'eff = 3.0: (dl) shows that after 5.2 mm a compressed few-cycle soliton 
forms (10 fs FWHM, sub-2 cycles), after which a Cherenkov wave (marked #1) is formed. It 
soon detaches from the soliton as its group velocity is much slower In the spectrum (d2) this 
Cherenkov wave has a peak at Adw = 2.65 /im. In Fig. [TJb) we use a longpass edge filter at 
2.5 jlm to isolate the Cherenkov wave: it is less than four optical cycles (35 fs FWHM) when 
formed, it is broadband (245 nm FWHM, or 350 cm^') and is near the transform limit. It then 
spreads out temporally but importantly the maximum peak intensity grows because the soliton 
keeps feeding radiation into the Cherenkov wave. This causes a GVM-induced broadening of 
the Cherenkov wave, which is in addition to that caused by pure GVD. The soliton then relaxes 
(decompresses) so at z = 8 mm it is too weak to feed radiation into the Cherenkov wave. At this 
point the Cherenkov wave detaches temporally from the soliton, and the observed drop in inten- 
sity is explained by dispersive broadening. Note that unlike the self-focusing case, where Ra- 
man red-shifting causes the soliton to slow down and continuously collide with the Cherenkov 
wave 1341 . here the Raman red-shift actually speeds up the soliton and the Cherenkov wave de- 
taches for good without any further interaction |14|. Figure[TJc) shows the filtered Cherenkov 
wave at different stages: it has and excellent near-Gaussian pulse quality, and the conversion 
efficiency starts at 4% and increases to over 10% with further propagation. Eventually soli- 
ton fission occurs 13511 creating a delayed minor soliton, which atz = 10—12 mm emits a 
new Cherenkov wave (marked #2) with a similar bandwidth and center wavelength as #1. This 
causes the mid-IR spectrum to split in two. The Raman-induced frequency red-shift of the major 
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Fig. 2. Spectral contents at the soliton formation point and at later stages on a log-scale 
for (a) Ai = 1.3 /im, (b) Ai = 1.65 flm and (c) A] = 1.8 ^m. In all three cases the initial 
pulse has 50 fs FWHM and A'eff = 3.0 was kept constant implying an input intensity of (a) 
/in = 215 GW/cm^ (b) Ii„ = 77 GW/cm^, and (c) Ii„ = 35 GW/cm^. 

soliton is larger than that of the minor soliton because its duration is shorter, so its Cherenkov 
resonance point should be closer to Xzd than that of the minor soUton. This is not the case, and 
we believe that it is because the spectral slope of the minor soliton is much steeper, pushing the 
Cherenkov wave to shorter wavelengths. At z = 15 mm the main soliton has recompressed and 
it emits a new burst of Cherenkov radiation, marked #3. An XFROG spectrogram lfT4l at z = 20 
mm is shown in (e), elucidating the temporal and spectral location of the interacting waves. 

In Fig. 12] the spectral evolution from Fig. [T] is compared to other pump wavelengths. In all 
cases the peak Cherenkov intensity grows upon propagation, and later several peaks appear be- 
cause new Cherenkov waves are emitted by soliton recompression and fission. When pumping 
far from Xzd, the Cherenkov radiation is very weak, see Fig. |2a) where Ai = 1.3 /xm, but in 
return it has 2-cycle duration. As the pump wavelength is increased, the Cherenkov wave shifts 
to lower wavelengths since the resonance wavelength changes, and its bandwidth reduces. As 
less bandwidth is required to sustain shorter-wavelength few-cycle pulses, a 4-cycle near-IR 
Cherenkov wave with Adw = 2.1 /xm can still be isolated in case (c) where Ai = 1.8 flm. The 
Cherenkov wave also becomes more intense when the pump wavelength approaches Azd be- 
cause the spectral overlap with the soliton is larger lfT9l . In (c) the pump is quite close to Xzd, 
and the red-shifted SPM shoulder leaks into the anomalous dispersion region (see also Ref. 
lfT4l ). This complicates the propagation dynamics after the initial soliton formation in a similar 
way as soliton fission. Despite this, the efficiency is very high, up to 25%. Finally, note that a 
longer interaction length is needed close to Azd because the dispersion length increases. 

4. Conclusion 

Trough numerical simulations of ultrafast noncritical cascaded SHG in LN we showed that 
few-cycle solitons can be formed that shed near- to mid-IR optical Cherenkov radiation in the 
A = 2.2—4.5 jim range with few-cycle duration, excellent pulse quality, and a high conversion 
efficiency (up to 25% was observed). We recommend a low soliton order to keep perturbations 
to the main soliton at a minimum so it can boost the Cherenkov wave over an extended prop- 
agation stage without distorting it. The Cherenkov waves could be isolated using a longpass 
filter, and 2-cycle mid-IR pulses were found when pumping far from the zero-dispersion wave- 
length Azd- This case has a low conversion efficiency, but it increases when pumping close to 
Azd where the Cherenkov wave shifts to shorter wavelengths, and the pulse duration increase 
slightly to 3-4 cycles when formed. Alternative methods for generating energetic few-cycle 
mid-IR pulses - like optical rectification ll36l or noncollinear optical parametric amplification 
ll37l - are either inefficient or very complex. Thus, cascaded SHG might provide an efficient 
bridge between near-IR femtosecond laser technology and ultrashort energetic mid-IR pulses. 



